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Abstract. We use orthogonality measures of Askey- Wilson polynomials to 
construct Markov processes with linear regressions and quadratic conditional 
variances. Askey- Wilson polynomials are orthogonal martingale polynomials 
for these processes. 



1. Introduction 

Orthogonal martingale polynomials for stochastic processes have been studied 
by a number of authors, see[51[Il[ISl[ll[Ml[151[lIl[lE]- Orthogonal martingale 
polynomials play also a prominent role in non-commutative probability [TJ , and 
can serve as a connection to the so called "classical versions" of non-commutative 
processes. On the other hand, classical versions may exist without polynomial 
martingale structure, see [6]. In [8] we identify intrinsic properties of the first two 
conditional moments of a process which guarantee that a stochastic process has or- 
thogonal martingale polynomials. These properties, linear conditional expectations 
and quadratic conditional variances, which we call the quadratic harness proper- 
ties, already lead to a number of new examples of Markov processes jHllinilll] with 
orthogonal martingale polynomials. Random fields with harness property were in- 
troduced by Hammersley [14] and their properties were studied e. g. in jlQ ) I32 j . 

In this paper we use measures of orthogonality of Askey- Wilson polynomials to 
construct a large class of Markov processes with quadratic harness property that 
includes previous examples either as special cases or as "boundary cases". The 
main step is the construction of an auxiliary Markov process which has Askey- 
Wilson polynomials ,4 as orthogonal martingale polynomials. The question of 
probabilistic interpretation of Askey- Wilson polynomials was raised in 12, page 
197]. 

The paper is organized as follows. In the reminder of this section we recall 
background material on quadratic harness property, Askey- Wilson polynomials, we 
state our two main results, and conclude with some applications. In Section [2] we 
give an elementary construction that does not cover the entire range of parame- 
ters, but it is explicit and does not rely on orthogonal polynomials. The general 
construction appears in Section [S] this proof follows the method from [9J and relies 
on martingale property of Askey- Wilson polynomials which extends projection for- 
mula for Askey- Wilson polynomials [20] to a large range of parameters. Section |4] 
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contains another elementary but computationally more cumbersome construction 
of a purely discrete quadratic harness which is not covered by Theorem ll.il 

1.1. Quadratic harnesses. In iS] the authors consider square-integrable stochas- 
tic processes on [0, oo) such that for all t,s > 0, 

(1.1) EiXt) = 0, E{XtX,) = mm{t, s}, 

E{Xt\J^s,u) is a linear function of Xs,Xu, and Var[Xf |J^s,ti] is a quadratic function 
of XgjXu- Here, J^s.u is the two-sided (T-field generated by {Xr : r e (0, s)U(w, oo)}. 
We will also use the one sided u-fields generated by {Xt r < t}. 
Then for all s < t < u, 

(1.2) E{Xt\T,.u) = ^X, + 

u — s u — s 

and under certain technical assumptions, [3 Theorem 2.2] asserts that there exist 
numerical constants rj, 9, a, r, 7 such that for all s < t < u. 



(1.3) Var[Xt|J-,, 



{u-t){t-s) , {uX,-sXuf uX, - sX^ 

^ I.. „\9 ' ^ 



u(l + as) -I- T — 7s \ (u- 
, i^u — XsY , ^Xu — Xs (Xu — Xs){uXs — sXj^ 



(u — s)2 ' u — s ^ — s)2 

In this paper it will be convenient to extend [8] by admitting also processes 
defined on a time domain T which is a proper subset of [0,oo), and by relaxing 
slightly condition (jl.ip . (Compare [SI Proposition 2.1].) Accordingly, we will say 
that a square-integrable stochastic process {Xt)t£T is a quadratic harness on T with 
parameters {r],9,a,T,^), if it satisfies (|1.2p and (|1.3I) on T. 

Under the conditions listed in [8l Theorem 2.4 and Theorem 4.1], quadratic 
harnesses on (0, 00) have orthogonal martingale polynomials. Although several 
explicit three step recurrences have been worked out in [3 Section 4] and for some 
of them corresponding quadratic harnesses were constructed in a series of papers 
[HI Uni E] J the general orthogonal martingale polynomials have not been identified, 
and the question of existence of corresponding quadratic harnesses was left open. 

It has been noted that the family of all quadratic harnesses on (0,oo) that sat- 
isfy condition (|l.ip is invariant under the action of translations and reflections of 
M: translation by a e R acts as (Xt) 1— e~'^Xf.-^af and the reflection at acts 
as (Xt) H> {tXi/t). Since translations and reflection generate also the symmetry 
group of the Askey- Wilson polynomials [21] , it is natural to investigate how to re- 
late the measures of orthogonality of the Askey- Wilson polynomials to quadratic 
harnesses. The goal of this paper is to explore this idea and significantly enlarge 
the class of available examples. We show that quadratic harnesses exist and are 
Markov processes for a wide range of parameters (77, 9, a, r, 7). The basic Markov 
process we construct has Askey- Wilson polynomials as orthogonal martingale poly- 
nomials. The final quadratic harness is then obtained by appropriate scaling and 
deterministic change of time. 

Since Askey- Wilson polynomials together with their limits are regarded as the 
most general classical polynomials, it is natural to inquire whether all quadratic 
harnesses can be obtained by taking limits of processes with Askey- Wilson transi- 
tion probabilities. 
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Theorem 1.1. Fix parameters —1 < q <1 and A,B,C,D that are either real or 
{A, B) or {C,D) are complex conjugate pairs, and such that ABCD,qABCD < 1. 
Assume that 

(1.4) AC, AD, BC, BD, qAC, qAD, qBD, qBD S C \ [1, oo). 

Let 



(1.5) r, 

(1.6) 

(1.7) a 

(1.8) r 

(1.9) 7 



{{A + B){1 + ABCD) - 2AB{C + D))^/l~^ 
^{l - AC){1 - BC){1 - AD){1 - BD){1 - qABCD) ' 

{{D + C)(l + ABCD) - 2CD(A + B))yT~~^ 
v/(l - AC){1 - BC){1 - AD){1 - BD){1 - qABCD) ' 
ABil - q) 



1 - qABCD ' 
CDjl - q) 
1 - qABCD ' 
q - ABCD 



1 - qABCD 

Consider the interval J = (To,Ti) of positive length defined by 



mm T L 7-1 + v/(7-l)^-4ar 
(1-10) To = max-^O, — , 



(1.11) — = max<^0, — ,-<J> ■ 

Then there exists a bounded Markov process (Xt)t(£j with mean and covariance 
(jl.ip such that p.2|) holds, and (|1.3p holds with parameters r],9,a,T,j. Process 
{Xt)ti£j is unique among the processes with infinitely- supported one- dimensional 
distributions that have moments of all orders and satisfy (|l.ll) . (jl.2p . and (|1.3p 
with the same parameters rj, 9, a, t, 7. 

The assumptions on A, B, C, D are dictated by the desire to hmit the number of 
cases in the proof but do not exhaust aU possibihties where quadratic harness {Xt) 
with Askey- Wilson transition probabihties exists; see Proposition 23] 

Remark 1.1. In Theorem ll.il (7 — 1)^ > 4aT. To see this, write (7 — 1)^ — 4(tt = 
{l~ABCD)'^{l-q)'^/{l-qABCD)'^. Together with the upper bound in [8, Theorem 
2.2] this shows that to construct a quadratic harness on [0, 00) from Theorem 1 1.1) 
we must have — 1 < 7 < 1 — 2y/aT. 

Remark 1.2. In terms of the original parameters, the end-points of the interval are: 

(1.12) To ^ niax{0,-CA^^^}, 

(1.13) 1 = max|o.-AS,^i^il^i. 
^ ' Ti \ ' ' 1 - qABCD J 

If CD < or AB < then the third term under the maximum contributes for 
q <0 only. 
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Remark 1.3. From (|1.9|) . the construction will give 7 > 1 when ABCD < — 1. 
Multiplying (|1.7p and ()1.8|) . we see that this may occur only when ctt < 0, i.e. for 
harnesses on finite intervals; compare [SI Theorem 2.2]. 

1.2. Martingale property of Askey- Wilson polynomials. For a,b,c,d S C 
such that a 7^ 0, and 

(1.14) abed, qabed, ab, qab, ae, qac, ad, gad G C \ [1, 00), 

the Askey- Wilson [4] polynomials Wn{x) = ijUn{x; a, b, c, d) are defined by the recur- 
rence 

(1.15) {2x - {a + a~^))wn{x) = A„iu„+i{x) - {A„ + C„)w„ -I- CnWn^i{x), n > 0, 
where 

(1 - afecdg"-i)(l - a6g")(l - acg")(l - adg") 



(1.16) An = 

(1.17) Cn = 



ail - abedq^''-^){l - abcdq^'') 
a(l - g")(l - 6cg"-i)(l - 6dg"-i)(l - cdg""!) 



(1 - afecdg2n-2)(i _ a5cdg2«-i) 

The initial conditions are 1 = and iDq = 1 . We also note that recurrence (|1.15p 
is symmetric in parameters 6, c, d but a is distinguished. 

Askey- Wilson polynomials 'Wn{x;a,b,c,d) can be defined also for a = 0; the 
definition appropriate for our purposes uses the limit as a —> of the normalized 
(monic) form of the three step recurrence: 

(1.18) XWn{x) = Wn+i{x) - i(A„ -f- C„ - (a + a^'^))Wn + ^An-lCnWn~l{x). 

Recall that polynomials {r„{x;t) : n G Z_|_,f G /} are orthogonal martingale 
polynomials for process {Zt)tei if^ 

(i) E [r,,{Zt;t)r„,{Zt;t)) = for m 7^ n and t G /. 

(ii) E \rnlZt;t)\Fs) = r'„(Zs; s) for s <tml and all n = 0, 1, 2 ... . 

The following result shows that Askey- Wilson polynomials are orthogonal mar- 
tingale polynomials for a family of Markov processes. If A, B, C, D are as in Propo- 
sition [^Hl this is a re- interpretation of the projection formula from pOj . 

Theorem 1.2. Suppose that A,B,C,D satisfy the assumptions in Theorem 
Let 

(1.19) / = I{A, B, C, D, q) = f max{0, CD, qCD}, ] 

\ maxjO, AB, qAB\ ^ 

with the convention 1/0 — 00. (The last terms under the maxima can contribute 
only when q < and CD or AB are negative.) If A ^ 0, let 

(1.20) r„(x; t) = u;„ {^ ^^"q'^ ^Vt, B^t, C/Vt, D/Vi^ . 

(If A = 0, use monic Wn instead ofwn-) Then 

{rn{x;t) : n = 0,1,2..., t G /} 

are orthogonal martingale polynomials for a Markov process (Zt) which satisfies 
()1.2p and ()1.3|) with r] = 9 = a = T = Q and j = q. 



Remark 1.4. Process {Zt) from Theorem ll.2l is closely related to the Markov process 
{Xt) that appears in the statement of Theorem I l.l) see (|2.27p . 
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Remark 1.5. If yl, B 7^ we can choose either A or B as the distmguished first 
parameter in ui„ . The corresponding polynomials r„ differ by the factor 'b^'iac'^d)" 
which does not depend on t, so both sets of polynomials are orthogonal martingale 
polynomials for the same process (Zt). 

1.3. Some worked out examples. This section shows how Theorem ll.ll is related 
to some previous constructions, and how it yields new examples. From examples 
that have been previously worked out in detail one can see that the boundary of 
the range of parameters is not covered by Theorem 1 1.1) in particular it does not 
cover at all the family of five Meixner Levy processes characterized by quadratic 
harness property in |31j . On the other hand, sometimes new examples arise when 
processes run only on a subinterval of (0, c»). 

1.3.1. q-Meixner processes. Theorem 11.11 allows us to extend |10i Theorem 3.5 ] 
to negative r. (The cases 7 = ±1 which are included in [10] are not covered by 
Theorem O) 



Corollary 1.3. Fix t, € M and — 1 < 7 < 1 and let 

(0 ifT>0 
ro=<-T/(l-7) «/t<0, 7>0 
{ T i/r < 0, 7 < 0. 

Then there exists a Markov process (Xt) on [To,oo) such that (jl.ip . (II. 2|) hold, and 



()1.3p holds with parameters rj — 0, a — 0. 
Proof. Let g = 7, A = 0, B = 0, and 



C = i — n = J 2Vl -<; — 

K 2^1-9 k 2^1-9 

Then p.4p holds trivially so by Theorem [TTTl {Xt) is well defined on J = [Tq, 00). 
Straightforward calculation of the parameters from (|1.6p , (jl.Sp , l|1.9p ends the proof. 

□ 

When r < 0, the univariate laws of Xt are the "sixth" family to be added to the 
five cases from pIJl Theorem 3.5 ]. The polynomials orthogonal with respect to the 
law of Xt satisfy recurrence 

xpn{x;t) ^pn+i{x;t) +0[n]qPn{x;t) + {t + T[n~ l]g)[n]qPn-i{x;t), 

where [n]q = {1 — q^^)/ {1 — q). So the polynomials orthogonal with respect to the 
standardized law of Xt/Vt are 

9 T 
(1.21) xpn{x;t) = pn+i{x;t) + —[n]qpnix;t) + (1 + -[n - l]q)[n]qPn-iix;t). 

The same law appears under the name g-Binomial law in [23) for parameters 
n = -t/r € N, r = -p{l - p) & [-1/4,0). When g < 0, this law is discrete 
two-point law at t = |r|. 

A justification of relating this law to the Binomial can be given for q = 0. In this 
case, recurrence (|1.2ip appears in 7, (3)] with their a = ^ and their b — j. By [3 

Proposition 2.1], the law I't of ':^^t is a free convolution -py-fold power of the two- 
point discrete law that corresponds to t = — t. That is, vt — I'-l'^^^', in particular. 
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at t = — nr, Xt/\/T has the law that is the n-fold free additive convolution of a 
centered and standardized two-point law. 

1.3.2. Bi-Poisson processes. Next we deduce a version of Theorem 1.2]. Here 
we again have to exclude the boundary cases 7 = ±1 as well as the case 1 + 776* = 
max{7, 0}. 

Corollary 1.4. For —1 < 7 < 1, and 1 + rjB > niax{7, 0} there exists a Markov 
process (-'''t)tiE[o,oo) such that (|l.ip . (jl.2p hold, and (jl.3l) holds with a = t = 0. 

Proof. LctA^O,B = -^^,C = 0,D = -^^^. Then ED = < 
1. The condition qBD < 1 is also satisfied, as we assume rjd + 1 > when q < 0. 
Thus (|1.4[) holds and we can apply Theorem ll.il From formulas (|1.5p through (|1.9I) , 
the quadratic harness has parameters ri,6,a — 0,t — 0, 7, as claimed. □ 

1.3.3. Free harness. Next we indicate the range of parameters that guarantee ex- 
istence of the processes described in [8l Proposition 4.3]. Let 

(1.22) a = 1±^,P^^. 

1 — err 1 — (TT 

Corollary 1.5. For < ar < 1, j — —ar, and r],d with 2 + rjO + 2(jt > and 
1 + a/3 > 0, there exists a Markov process (^t)te[o.oo) such that (jl.ip . (jl.2p and 
(1131) hold. 

Remark 1.6. When 2 + rjB + 2aT < 0, two of the products in (|1.4p are in the 
"forbidden region" [l,oo), so Theorem 1 1 . II does not apply. However, the univariate 
Askey- Wilson distributions are still well defined. 

Proof of Corollary \1.5[ Take q — 0, and let 

+ I3a- \J~A cr+{a- /3crf a + /3 a + \J-Aa+ {a - P af 



A = ' , , B 



2^1 + af3 ' 2Vl + afi 



(3 + aT- ^-4:T+ {P-arf jS + ar + y^-4 t+{I3- arf 

~ 2Vl + a/3 ' ^ ^ 2 VI + ' 

To verify that AC ^ [l,oo) we proceed as follows. Note that 

a + (t/3 



(1.23) A + B = 

(1.24) C-hZ? = 

(1.25) A-B = 

(1.26) C~D ^ 



VI + a/5 
ar + /3 

" VI + a/3 ' 
^(a - cr/^)2 - 4ct 

VI + a,3 
^(/3 - ra)2 - 4t 



Vl + a/i 

Multiplying (A -f- B)(C -1- D) and (A - B){C - D) and using ABCZ? = err we get 



^^^AC ^^^BC^ TT^ 

and 

AC BC 1 + aj3 



ASKEY-WILSON POLYNOMIALS, QUADRATIC HARNESSES AND MARTINGALES 7 



This gives the foUowing quadratic equation for AC: 



(1.27) + — - 2(1 + a/3) • 

We now note that a quadratic equation x + aj x — h with < a < 1 and complex 
h can have a root in [1, oo) only when h is real and 6 > 1 + a; this follows from the 
fact that X + ajx is increasing for a; > a, so x + a/x > 1 + a for a; > 1. 

Suppose therefore that the right hand side of (|1.27p is real and larger that l + crr. 
Then calculations lead to — Aa^/O'^ — 4t > 2 + r]6 + 2aT. The right hand side is 
nonnegative by assumption, so squaring the inequality we get (1 + a /3) (1 — cr t) < 
0, which contradicts the assumption. 

Other cases with AD, BC, BD are handled similarly. Since ABCD = ctt < 1 by 
assumption, by Theorem 1 1 . 1 1 the quadratic harness exists. 

It remains to calculate the parameters. From AB = a CD = t we see that 
()1.7|) and ()1.8p give the correct values, and 7 = — ctt from ()1.9|) . To compute 
the remaining parameters, we re-write the expression under the square root in the 
denominator of (11.51) as 



(1 - AC)il - BC)il - AD){1 - BD) 

^(l + .r-(AC+^)) {l + <rr-iBC+^) 

This is the product of two conjugate expressions (see (|1.27p . and its derivation). A 
calculation now simplifies the denominator of (|1.5|) to (1 — ar) / \/T~\~a^ . Inserting 
([L231) and (fr24| . the numerator of (fTSl) simplifies to {a - l3a) (1 - (tt)/V1 + a/3 
The quotient of these two expressions is a — f3a — rj. Similar calculation verifies 

(iTTel). □ 



1.3.4. Purely quadratic harness. The quadratic harness with parameters 77 = = 0, 
and (TT > has not been previously constructed. 

Corollary 1.6. For cr, r > with ctt < 1 and —1 < 7 < 1 — there exists 

a Markov process (^t)tg[o,oo) such that (|1.1|) . ()1.2|) hold, and (|1.3|) holds with rj = 

Proof. Let 

4(7 + err) 



1 + 7 + y (1 - 7)^ - 4crr 
To see that — 1 < g < 1, note that for 7 + ctt ^ 0, 



1 +7^ - 2ctt - (1 + 7) y (1 - 7) - 4crT 
2 (7 + aT) ' 



which gives 



-2a/(1 -7)^ -4crr 



(1.28) q-l = , = < 



l + 7+\/(l-7)^-4crT 
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and 

(1.29) q + l = ' /^ + ^) 

1+7+ \/(l- 7)^-4 CTT 

Noting that (1 - q)^ + iqar > 4ctt(1 - ar) > 0, let 



^{l-q) + ^{l-qf +4qa' 



and 

C = -D 



\j{l-q) + ^{l-qf +AqnT 

Since B, C, D are purely imaginary, we only need to verify condition BC < 1 
which reads 



(1.30) q + 2y/^ -l<^J{l-qf+ 4qaT. 

This is trivially true when q + l^far — 1 < 0. If g + l^/ar — 1 > 0, squaring both 
sides we get 4(1 — q)y/aT > 4(1 — q)a-T, which holds true as q > 1 and < err < 1. 

Thus quadratic harness {Xt) exists by Theorem 1 1.1) and it remains to verify that 
its parameters are as claimed. A straightforward calculation shows that (|1.7p and 
(|1.8[) give the correct values of parameters. It remains to verify that formula (|1.9p 
indeed gives the correct value of parameter 7. Since this calculation is lengthy, we 
indicate major steps: we write (|TT9|) as7-l = {q - + ABC D) /{I - q ABC D), 
and evaluate the right hand side. Substituting values of A, B, C, D wc get 



(, - 1)(1 + ABCD)/il - qABCD) = - + + ^<I-r _ 

2q 

Then we use formulas (|1.28p and (I1.29P to replace 1 — q and I + q and note that 
since 7 < 1 — 2y^aT we have 7 < 1 — 2(tt and 

r, 2(1 - 7 - 2(tt) 

V(l-g)2+4gaT- \ \ ' =. 

7+ V(l - 7) - 4(TT+ 1 

This eventually simplifies {q- ABC D) /{I- qABCD) to 7-I, as claimed. □ 

2. The case of densities 

In this section we give an explicit and elementary construction of a quadratic 
harness on a possibly restricted time interval and under additional restrictions on 
parameters A, B, C, D. 

Proposition 2.1. Fix parameters —1 < g < 1 and A,B,C,D that are either real 
or (A, B) or (C, D) are complex conjugate pairs. Without loss of generality, we 
assume that \A\ < \B\, \C\ < \D\; additionally, we assume that \BD\ < 1. Then 
the interval 

J \\D\'^-CD I-CDIBI"^' 



1-AB\D\^' \B\^-AB 
has positive length and there exists a unique bounded Markov process {Xt)teJ with 
absolutely continuous finite dimensional distributions which satisfies the conclusion 
of Theorem ] 1. 11 
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Remark 2.1. Proposition 12 . 1 1 is a special case of Theorem I l.lj the latter may allow 
to extend the processes constructed here to a wider time-interval: 

\D\^-CD , 1-CD\B\^ 

The easiest way to see the inequalities is to compare the end points of the intervals 
(f2^ and ([LTQII after the Mobius transformation (|2?26l) : for example, > \CD\ > 
max{0, CD, qCD}, where the last term plays a role only when q < and CD < 0. 

The rest of this section contains the construction, ending with the proof of Propo- 
sition 12.11 



" ^ ^ aqi) rt = 1, 2, 



2.1. Askey-Wilson densities. For complex a and \q\ < 1 we define 

{a)n = {a\q)n = 




n = Q 



(a)oo = (a; q)oo = n (1 - 

and we denote 



(ai,a2, . . . ,a;)oo = (01,02, . . . ,a/;g)oo = (oi; 9)00(02; g)oo • ■ • {ai;q)oc , 

(ai, 02, . . . ,a;)„ = (oi, a2, . . . , o/; g)„ = (oi; 5)00(02; g)oo ■ • • {ai;q)n ■ 
The advantage of this notation over the standard product notation lies both in its 
conciseness and in mnemotechnic simplification rules: 

(a, b)n _ {b)„ 
(a, c)„ (c)„ 

(2.1) ia)M+L = (9*^a)L (a)M , 
and 

(2.2) {a)M = {~af\^^ ' ^ 



which often help with calculations. For a reader who is as uncomfortable with this 
notation as we were at the beginning of this project, we suggest to re- write the 
formulas for the case q = 0. For example, (o;0)„ is either 1 or (1 — a) as n = or 
n > respectively. The construction of Markov process for g = in itself is quite 
interesting as the resulting laws are related to the laws that arise in Voiculescu's 
free probability; the formulas simplify enough so that the integrals can be computed 
by elementary means, for example by residua. 

From Askey and Wilson [4, Theorem 2.1] it follows that if a,b,c,d are complex 
such that max{|a|, |&|, |c|, \d\} < 1 and — 1 < q < 1, then with 9 — 6^ such that 
cos 9 — X 

(•2 3) / -^=^= ■ i-2£ dx 

2i:(abcd)aa 



{q, ah, ac, ad, be, bd, cd)c 
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When —l<q<l and a, b, c, d are either real or complex and form two conjugate 
pahs and max{|a|, |c|, \d\} < 1, the integrand is real and positive. This allows 
us to define the Askey- Wilson density 

2 



(2.4) /(x;a,6,c,d) 



K{a, b, c, d) 



where 

(q, a6, ac, ad, be, bd, cd)o 
(2.5) -fi (a, 0, c, a) = — 



A-i.i)(^) ' 



2n{abcd)oo 

The first two moments are easily computed. 

Proposition 2.2. Suppose X has distribution f(x; a, b, c, d) with parameters a, b, c, d 
as above. Then the expectation of X is 

, , wf/v s a + b + c + d ~ abc — abd — acd — bed 

(^•^^ ^(^^ = 2{l- abed) ' 

and the variance of X is 

(1 - ab){l - ae){l - ad){l - 6c) (1 - bd){l - ed){l - q) 



(2.7) Var(X) = 



4:{l-abed)^{l-abedq) 



Proof. lfa = b = e = d = 0, E{X) = by symmetry. If one of the parameters, say 
a G C is non-zero, we note that 



/ ie -ie\ I ie -ie 
[ae , ae ' — / 



\ f ie -ie\ ( ie -ie\ 
)oo = («e ' «e )i He ,aqe )^ 

= {l + a'^~2ax) {aqe'^,aqe~'^) 



2 n V, K{a,b,e,d) {I - ab){l - ae){l ~ ad) 



Therefore by ([M 

E(l + a^ - 2aX, , , n 

K[qa, 0, e,d) 1 — abed 

Now (|2.6p follows by a simple algebra. 
Similarly, for non-zero a, 6 G C, 

4a6 Var(X) = E [(1 + - 2aX)(l + 6^ - 26X)] -E(l + a2-2aX)E(l + 62-26X) 

K{a,b,c,d) K^{a,b,e,d) 
K{qa, qb, c, d) K{qa, b, c, d)K{a, qb, e, d) 
_ (1 - a6)(l - qab){l - ae){l - ad){l - be){l - bd) 
(1 — abed){l — qabed) 

(1 - a6)2(l - ac)(l - arf)(l - 6c) (1 - bd) 
(1 - abed)^ ' 

Again after simple transformations we arrive at (j2.7p . 

If only one parameter is non-zero but q ^ 0, the calculations are similar, start- 
ing with E ((1 + a'^ — 2aX){l + a^q^ — 2aqX)); when q = the density is a re- 
parametrization of Marchenko-Pastur law [121 (3.3.2)]; we omit the details. If 
a,b,c,d are zero, /(x; 0, 0, 0, 0) is the orthogonality measure of the continuous q- 
Hermite polynomials [T71 (3.26.3)]; since H2{x) = 2xHi{x) - (1 - q)Ho = ^x^ - 
[\ — q), the second moment is (1 — (7)/4. □ 

We need a technical result on Askey- Wilson densities inspired by (301 formula 
(2.4)]. 
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Proposition 2.3. Let a, b, c, d, q be as above with the additional assumption that 
the only conjugate pairs possible are a — b and/or c — d and m is real such that 
\m\ < 1. Then with x = cos Ox, 

(2.8) J f {x;am,bm,c,d) f (^y^ajb^me^^" ,me~^^'') dx — f {y;a,b,cm,dm) . 

Proof. We compute the left-hand side of (|2.8p expanding the constants K{am, bm, c, d) 
and A'(a, 6, me*^"^ , me~*^"^) to better show how some factors cancel out. To avoid 
case-by-case reasoning when complex conjugate pairs are present, we also expand 
parts of the density without the use of modulus as in (|2.3p . 
The integrand on the left hand side of (|2.8|) is 



(q, abm^ , acm, adm, bcm, bdm, cd)ao \ i^'^^^")^ P 
2Tr (abcdm"^) oo (ame^^^ , ame~**^, bme'^^^ , 6me^*^"=, )^ | (ce*^»', de*^"^)^^ P 

(g, ab, ame'^=", ame~*^=", &me*^^, &me~*^^, m^) 
X , = 



Rearranging the terms we rewrite the left hand side of p.Sp as 



(q, abm? , acm, adm, bcm, bdm, cd, q, ab, m^)oo I (^^*^'')oo 1^ 



X 



(2^)2(a6m2, abcdm^)^./T^ I (ae*''", be'^y)^ P 

|(e-«^)^P dx 



Now we apply formula (12. 3p to this integral, so the left hand side of (j2.8p becomes 
(g, abm^, acm, adm, bcm, bdm, cd, q, ab, m'^)oo I (e^*^""* 



X 



27r(cdm2)o£j 

{q, m2, mce*^", mde^^y, mce^^^y , mde^^^y, cd)co 



{q, ab, acm, adm, bcm, bdm, cdm'^)oo (^^*^'')oo 1^ 

X 



which completes the proof. □ 

2.2. Markov processes with Askey- Wilson densities. We now fix A, B, C, D 

as in Proposition 12. II The interval 

1 

is non-empty (here 1/0 = oo). For any t E I{A, B, C, D) and y E [—1, 1] let 
(2.10) p{t,y) = f(^y;AVt,BVi,-^^,^^ 



(2.9) IiA,B,C,D)= [\D 
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and for any s < t in I {A, B, C, D) and x,y £ [—1,1] let 

(2.11) p(s,x;^,y) = /(^2/;AV^,B^/^,^e'^^^e-'^==^ , x = cos 0, . 

Proposition 2.4. The family of probability densities {p{s,x;t,y), p{t,y)) defines 
a Markov process (i^t)tg/ on the state space [—1, 1]. That is, for any s < t from 
I{A,B,C,D) andy e [-1,1], 

(2.12) p(t,y) ^ J p{s,x;t,y)p{s,x) dx 
and for any s < t < u from I{A, B, C, D) and x, z £ [—1, 1], 

(2.13) p{s,x;u,z)= p{t,y;u,z)p{s,x;t,y) dy . 



Proof. To show (|2.12p it suffices just to use the identity (|2.8p with a = Ay/t, 
b — B\/t, c — C/y/s, d ~ D/^/s, and m — \fsjt G (0,1). We note that this 
substitution preserves the conjugate pairs in the case of complex parameters, and 
that by the definition of I{A, B,C, D), parameters Ay/i, By/t, C/y/s and D/y/s 
have modulus less than one, so ()2.8p applies here and gives the desired formula 



^f(y;AV-t,BV-t,^^,^^ 

To get the second formula (12. 13^ we again use (j2.8p this time with a — Ay/u, 
b — y/u, c = y/sjte^^=' , d = y/sjte~^^'' , and m — y/tju. Thus we arrive at 



f (z;AV^,BV^, ^e^«^ ^e-^««) / (^y; AVt, BVt, ^e'«^ ^e^^^) dy 



□ 



Proposition 2.5. Let {Yt)t,=i{A,B.c.D) be the Markov process from Proposition \2.4\ 
with marginal densities (|2.10p and transition densities (|2.1ip . Fort G I {A, B, C, D) 

(2 U) E(Y ) = + B - ABjC + D)]t + C + D - CDiA + B) 
^ ' ' ^ 2Vt(l-ABCD) 

(2.15) Var(y() 

_ (1 - g)(l - AC)il - AD)il - BC){1 BP) 

U{l-ABCDY{l^qABCD) ^^^^^ ' ^^'^ 
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and for s,t E I [A, B, C, D), such that s <t, 

(2.16) Cov(n,yt) 

(1 - q){l - AC){1 - AD)(1 - BC)(1 - ED) 



(2.17) E(rt|j-,) 

(2.18) Ya.r{Yt\Fs) 

(1 - q){t - s)(l - ABt) 



4v^(l - ABCDfil - qABCD) 

{A + B){t~s) + 2(1 - ABt)y/^Y, 



{s-CD){l- ABt) 



At{l - ABsf{l - qABs] 



2^ft{l - ABs) 



(1 + A^s - 2A^/:sYs){l + B^s - 2ByGY,) 



Proof. Formulas (|2.14l) and (|2.15p follow, respectively, from (|2.6p and (|2.7p by tak- 
ing a = AVi, b = By/i, c = C/^/^ and d = D/^/t. 

Similarly, the formulas (|2.17l) and (I2.18P follow, respectively, from (12. 6p and (I2.7| 
by taking a = A^ft, b = By/t, c = e'**- and d = e"'"- . 

To obtain the covariance we make use of ()2.17p as follows: 

Cov(n,yt) = EiYMYtlTs)) -^YsMYt) 

f{A + B)(t-s) , (l-ARs)yi 



V2^/^(l - ABs) 

Now the formula (|2.16p follows, after a little algebra, from (|2.14l) and (|2.15p . 



□ 



Next we show that the conditional distributions of Yt given the past and the 
future of the process is given by an Askey- Wilson density that does not depend on 
parameters A, B, C, D. 

Proposition 2.6. Let (Yt)t<=i(A,B.c.D) be the Markov process with marginal densi- 
ties ()2.10p and transition densities (|2.1ip . Then for any s, t, u G I{A, B, C, D) such 
that s < t < u the conditional distribution ofYt given Ts,u has the Askey-Wilson 
density 

(2.19) / (y; exp(i6'^), exp{-ie,), ^ exp(i6'a;), ^ exp(-i6'^) 

(Here, x = cos9x = Y^, z = cosOz = Y^-) The first two conditional moments have 
the form 

(2.20) E(r.i^..)^ ^"-^^-^;:f ^^r^^^" , 

\'t(U — s) 

(2.21) 

V'^YtlTs «) = - " -t)it-s) fl {u^Ys - s^Y^){^Yu -V^sYsY 



t{u ~ qs) \4 (u — s)^ 
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Proof. By the Markov property it follows that the conditional density is 

p{t,y;u,z)p{s,x;t,y) 
p{s,x;u,z) 

f{z; A^, B^, ^ e-^ , ^ e-"^« )/(?/; AVi, BVt, ^ e'^- , ^ g-'^- ) 
/(z; AV^, ^ e'S^ ^ e-^''^) 

Now the result follows by plugging in the formula above the definition of the Askey- 
Wilson density (|2.4I) with suitably chosen parameters. The mean and variance are 
calculated from (gH) and dUT]). □ 

Proof of Provosition \2. 11 If we define a new process {Zt)t£i(A,B.c.D) through 
(2.22) Z, = -^Y, , 

then is Markov and satisfies 

so that 



( AB,/T^Zt-{A + B) 
\ l-ABt ' * 



-(1 + A^s - A^\-qZs){l + B^s - B^l - qZs 



is a martingale. Moreover, 

Var(Zt|J-,) 

{t-s){l-ABt) 
~ (1 - ABs)2(l - qABq)' 

For the double conditioning with respect to the past and future jointly it follows 
that (Zt) satisfies quadratic harness conditions 

U — t t — 5 

(2.23) E{Zt\J^s.u) = Zs + Z„ 

u — s u — s 

and 

(2.24) Var(Z,|^,„) = (l~{l~q) ^"^^ '^^^^^^ " 



u — qs \ i"^ ^ s)'^ 

which correspond to the g-Brownian motion when the process is defined on (0, oo) 
and has covariance min{s,i}, see [lOl Theorem 4.1]. Here, (Zt) is defined only a 
finite time domain I{A, B, C, D) and the covariance is different: 

(2.25) Cov(Z., Z.) = i^-fC)il-AD)il-BC)il-BD) 

^ ' V t-- (I- ABCDY{1- qABCD) ^ ' 

(The law of Zt will differ from the g-Gaussian law if \A\ + \B\ + \C\ + \D\ > 0.) 

The covariance is adjusted by a suitable deterministic time change. Consider a 
Mobius transformation 

(2.26) M.) = ff^, 
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which for ABCD < 1 is an increasing function with the inverse T{t) = (t+CD) /(1 + 
ABt). Note that J = J{A,B,C,D) = h{I{A,B,C,D)). For t e J{A,B,C,D), 
define 

(2.27) Xt := Xt-A,B.c.D.q 

_ ZT(t) - ]E(^T(t)) ^ (1 - ABCD)^\ - qABCD 



1 - T{t)AB ^(1 - AC)il - SC)(1 - AD){1 - BD) 

_ VT^{\ + ABt)ZT(t)^{A + B)t-(C + D) ^^^^ 

v/(l - g)(l - AC){1 - AD){1 - BC){1 - BD) ^ 

A calculation shows that [Xt)t£j has unconditional and conditional moments as 
claimed: formula is a consequence of (|2.25p . and (|1.2p follows from (|2.23p . 

A much longer calculation shows that (|2.24p translates into (11.31) with parameters 

(HMm. □ 

3. Proof of Theorem 11.11 

Next, we tackle the issue of extending the quadratic harness from Proposition 
12.11 to a larger time interval. The main technical difficulty is that such processes 
may have a discrete component in their distributions. The construction is based 
on the Askey- Wilson distribution [4, (2.9)], with slight correction as in [29', (2.5)]. 

The basic plan of proof of Theorem 11.11 is the same as that of the proof of 
Proposition 12.11 we define auxiliary Markov process {Yt)tei through a family of 
Askey- Wilson distributions that satisfy the Chapman-Kolmogorov equations. Then 
we use formulas ()2.22p and (|2.27p to define [Xt). The main difference is that due 
to an overwhelming number of cases that arise with mixed type distributions, we 
use orthogonal polynomials to deduce all properties we need. (Similar approach 
was used in ^].) 

3.1. The Askey- Wilson law. The Askey- Wilson distribution iy{dx;a,b,c,d) is 
the orthogonality measure of the Askey- Wilson polynomials {u)„(a;; a, 6, c,d)}. Since 
{An} and {C'n} are bounded sequences, v{dx; a, b, c, d) is unique, and compactly 
supported [HI Theorem 2.5.4]. If |a|, |5|, |c|, \d\ < 1, this is an absolutely continuous 
measure with density (|2.4p . In general, v{dx; a, 5, c, d) may have a discrete compo- 
nent, or be purely discrete as in (|4.ip . We need also to point out that although 
parameter a is distinguished in the definition of w„, the orthogonality measure is 
symmetric in all parameters. This can be seen from representation 

(3.1) «;„(a:;a,6,c,d) = 4^3(^ 

and the 4(p3-Sears transformation [26] which shows how to rescale the Askey- Wilson 
polynomials to be independent of the order of parameters a, 6, c, d without affecting 
the orthogonality measure. 

In general, it is quite difficult to give explicit conditions on a, b, c, d for the 
existence of v{dx]a,b,c,d). By Favard's theorem [121 Theorem 2.5.2], measure 
v{dx; a, 6, c, d) exists for all a, 6, c, d such that Ak,Ck G R for all k and nfe=o "^kCk+i > 
for all n. Furthermore, it is known that if TV is the first integer such that 
An^iCn = 0, then u{dx;a,bjC,d;q) is a discrete measure supported on the finite 
set of N zeros of the polynomial wn{x). 
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Under (|1.14|) . if a, b, c, d are either real or come in complex conjugate pairs and 
a then Ak > and C/c £ R for all fc; so in this case the Favard condition is 

n 

(3.2) n '^'^ - ^ '"^^^ 

fc=i 

A simple sufficient condition for p.2p is that (Il.f4p holds and in addition 

(3.3) be, qbc, bd, qbd, cd, qcd € C \ [1, co) . 

Under this condition, if a, 6, c, d are either real or come in complex conjugate pairs, 
then the corresponding measure of orthogonality v{dx] a, b, c, d) exists. Unfortu- 
nately, this simple condition is not good enough for our purposes - we need to 
allow also Askey- Wilson laws with finite support as in 3 . In fact, such laws de- 
scribe transitions of the Markov chain in the atomic part. 

We now state conditions that cover all the cases needed in this paper. Let 
mi — mila, b, c, d) denote the number of the products ab, ac, ad, be, bd, cd that fall 
into [1, oo), and let m2 — m2{a, b, c, d) denote the number of the products qab, qac, 
qad, qbc, qbd, qcd that fall into [l,cx)). 

Lemma 3.1. Assume abed < 1 and qabed < 1. If q > then either mi — 
or mi — 2; if q < and mim2 = then there are only three possible ehoices: 
(mi, 771,2) = (0, 0), (0, 2), (2, 0). Only the following cases are possible: 

(i) If Q ^ and mi is zero, then u{dx; a, b, c, d) exists. 

(ii) If q < and mi = rn2 = 0, then u{dx; a, b, c, d) exists. 

(iii) If q > and mi = 2, then a,b,c,d S M. In this case, v[dx;a,b,c,d) is 
well defined if either q — or the smaller of the two products that fall into 
[1, 00) is of the form 1/q^ and in this latter case u{dx; a, b, c, d) is a purely 
discrete measure with N + \ atoms. 

(iv) If q < and mim2 > then v{dx; a, b, c, d) fails to exists. 

(v) If q < and mi = 2, m2 — then a, b,c,d € M. In this case, v{dx; a, b, c, d) 
is well defined if the smaller of the two products in [l,oo), say ab, is 1/q^ , 
with even N . Then v{dx; a, b, c, d) is a purely discrete measure with + 1 
atoms. 

(vi) If q < mi — 0, 7712 — 2, then a, b,c,d G M. In this case, v{dx; a, b, c, d) is 
well defined if the smaller of the two products in [l,oo), say qab, is 1/q^ , 
with even N . Then v[dx; a, b, c, d) is a purely discrete measure with N + 2 
atoms. 

Proof. The fact that ?Tii = or mi = 2 is an elementary observation. Existence 
of v{dx;a,b,c,d) follows by inspection of YYl^iAk-iCk- If m2 > then in cases 
(iii) and (v) 0^=1 ^fe-iC'fc > for n < A^ and is zero for n > Af -|- 1. In case (vi), 
HLi Ak-iCk = for ah n > A^ + 2. 

□ 

According to Askey and Wilson [4] the orthogonality law is 

iy{dx;a,b,c,d) = f{x;a,b,c,d)li^l<i + ^ p{x)Sx; 

xeF(a,b,c,d) 

here F = F{a, b, c, d) is a finite or empty set, and the density / is given by (|2.4|) . The 
continuous component is absent when K{a,b,c,d) = (recall ()2.5|) ). The atoms 
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occur for parameters which are real and outside of [—1, 1]. Each such parameter, 
say a, leads to a finite set of atoms at 

(3.4) = ^-sl±pn:l, 

with j >0 such that |<?-'a| > 1 and the corresponding probabilities are 

{a^'^,bc,bd,cd)ca 

3.5 p[xo) = — — 7—7—, 

(o/a, c/a, a/ a, abcdjaa 

(3 6) p{x ) = p{x^) ("^»^'°g'«'^)i(l ~ Q^g'^^) f 1 y" 

{q,qa/b,qa/c,qa/d)j{l — a^) \abcdJ ' ~ 

(This formula needs to be re- written in an equivalent form to cover the cases when 
abed = 0. It is convenient to count as "atom" the case \q^a\ = 1 even though the 
corresponding probability is 0.) 

We remark that if X has distribution i^{dx;a,b,c,d) then formulas for K{X) 
and Var(X) from Proposition 12 . 2 1 hold now for all admissible choices of parameters 
a, b, c, d, as these expressions can equivalently be derived from the fact that the first 
two Askey- Wilson polynomials integrate to zero. 

3.2. Construction of Markov process. Recall / — I {A, B, C, D; q) from (|1.19p . 
As in Section[2J we first construct the auxiliary Markov process {Yt)t^i. We request 
that the univariate law irt of Yt is the Askey- Wilson law 



(3.7) nt 



(dy) = 1^ (^dy; Ayft, BVt, 



In order to ensure that this univariate law exists, we use condition p.3p . This 
condition is fulfilled when (II. 4p holds and the admissible range of values of t is the 
interval / from (|1.19p . (The endpoints (|1.12l) and (11.131) were computed by applying 
Mobius transformation (|2.26p to the endpoints of /.) 

For t £ I, let Ut be the support of Trt{dy). Under the assumption (|1.4|) . this 
set can be described quite explicitly using the already mentioned results of Askey- 
Wilson [4]: Ut is the union of [—1, 1] and a finite or empty set Ft of points that are 
of the form 

(3.8) ^ 1 (bVI,' + ^) or = i + ;^ 

(3...) y,(t) = 1 [aV-,^ + o. MO = + 5?, 

There is at most finite number of points of each type. However, not all such atoms 
can occur simultaneously. All possibilities are listed in the following lemma. 

Lemma 3.2. Under the assumptions of Theorem without loss of generality 
assume \A\ < \B\ and \C'\ < \D\. Then the following atoms occur: 

• Atoms Uj{t) appear for D,C G M., and t E I that satisfy t < ; allowed 
j > satisfy D'^q^-' > t. 

• Atoms Vj{t) appear for Z),C G M, and t E I that satisfy t < ; allowed 
j >0 satisfy C'^q^^ > t. 

• Atoms Xj{t) appear for A, _B G K, and t E I that satisfy t > 1/ B^; allowed 
j > satisfy tB'^q^^ > 1. 
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• Atoms yj{t) appear for A, B G M, and t ^ I that satisfy t > l/A^ ; allowed 
j >0 satisfy tA^q^^ > 1. 

(The actual number of cases is much larger as in proofs one needs to consider all 
nine possible choices for the end points of the time interval /.) 
Next, we specify the transition probabilities of Yt. 

Proposition 3.3. For s <t, s,t €z I and any real x €z Us measures 

Ps,t{x, dy) = V (^dy; AVt, BVt, -^(x + Vx^), ^-{x - ^ x"^ - 1)^ 
are well defined. 

Proof. Clearly, Pg^t{x,dy) for a; G [—1,1] is absolutely continuous and has density 
(12.111) . This covers all possibilities when {A,B) and {C,D) are conjugate pairs or 
when l^iyi, |B|Vs, \C\/^/s, \D\/y/^ < 1, as then Us C [-1, 1]. 

It remains to consider x in the atomic part of TTs{dx). Relabeling the parameters 
if necessary, we may assume \A\ < \B\ and \C\ < \D\. For each of the cases 
listed in Lemma 13.21 we need to show that the choice of parameters a = A\Jt^ 



b = B^/t, c = y^{x + ^x^ — 1), d = ^J\{x ~ y/x^ — 1) leads to non- negative 
products rifc^o ^kCk+i > (recall (|1.16p and ([TTT71)). We check this by considering 
all possible cases for the endpoints of / and all admissible choices of x from the 
atoms of measure tt^. In the majority of these cases, condition p.3p holds so in fact 
AkCk+i > for all k. 

Here is one sample case that illustrates what kind of reasoning is involved in 
the "simpler cases" where p.Sp holds, and one example for more complicated case 
where ()3.3|) fails. 

• Case CD < 0, AB <0, q>0: In this case. A, B, C, D are real, / = (0, cxi) 
and assumption ABCD < 1 implies < 1/B^. A number of cases arises 
from Lemma 13.21 and we present only one of them: 

— Sub-Case x = Vj{s): Then < s < and the Askey- Wilson pa- 
rameters oi Ps^t{x]dy) are 

a — AVt, b = BVt, c — — d ~ 



Vt ' CqiVt 
Thus 

As Bs 
ab = ABt < < 1, ac = ACq\ ad = — — , be = q^BC, bd = — — , cd s/t < 1 

Cq^ Cq^ 

and 

qab = qABt <0<1, qac = ACq^ + ^ < ACq^ , qad = qbc = q^+^BC, 

Cq^~^ 

, , Bs 

^ Cqi-^ ' ^ ' 

Since s < C^q^^ , this implies \ad\ < \A\Vs < 1 as s < < < 
l/B"^ < \/A^. For the same reason, |6c?| < 1. Of course, \qad\ < \ad\ < 
1, \qbd\ < \bd\ < 1 

Finally, since AC,qAC, BC,qBC < 1, we have ac, qac, be, qbc < 1. 
Thus by Lemma l3.1f i). Ps^t{x,dy) is well defined. 
(We omit other elementary but lengthy sub-cases that lead to p.3p .) 
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• Case A,B,C,D eR,AB> 0, CD > 0: Rcyc I{A, B,C, D) = {CD,1/{AB)) 
is non-empty. Again a number of cases arises from Lemma [321 of which we 
present only one. 

- Sub-Case x = x.j{s): Here B'^ < s < t < 1/{AB) and B^\q^ > 1. 
Then the Askey- Wilson parameters of Ps.t{x; dy) are 



AVt, b = SVi, c = 



Here, Lemma l3.ir ii) or Lemma l3.ir v) applies with mi = 2 and TO2 = 
when 5 > or J is even, and Lemma I3.1f vi) applies with mi = 0, 
TO2 = 2 when g < and j is odd. To sec this, we look at the two lists 
of pairwise products: 



A 

ab = ABt < 1, ac = ABsq' < 1, ad = -, be = sq^B^, bd = 1/q^, 

Bqi 

cd = s/t < 1 

and 

A 

qab = qABt < < 1, qac = ABsq^^'^ < 1, qad = „ . , , qbc = sq^^'^B'^, 

Bq^ ^ 

qbd = l/g^ \ qcd = qs/t < 1 . 

Since |^| < \B\ we sec that l/iAB) < l/A^ so |^| < l/^/s and 

< < 1. This shows that ad < 1 and qad < 1. 

It is clear that both be, bd > 1 when q > or j is even, and that 
qbe, qbd > 1 when q < and j is odd. Thus by Lemma [Ol Pg^tix, dy) 
is well defined. 

Other cases are handled similarly and are omitted. □ 

In the continuous case, pt{dy) = p{t, y)dy and Pg^tix, dy) = p{s, x; t, y)dy corre- 
spond to (|2.10p and (|2.11[) respectively. We now extend Proposition 12 .41 to a larger 
set of measures. 

Proposition 3.4. The family of probability measures ntidy) together with the fam- 
ily of transition probabilities Ps,t{x, dy) defines a Markov process {Yt)tei on IJjgj Ut- 
That is, Ps^t{x, dy) is defined for all x £ Us- For any s < t from I{A, B, C, D) and 
Borel set F C M, 



(3.10) Tlt{V) = / Ps,t{x,V)Tis{dx) 

and for s < t < u from I, 



(3.11) PsAx. V) = / PtAy, V)Ps,t{x, dy) for all x e Us . 

JUt 

To prove this result we follow the same plan that we used in pj and we will 
use similar notation. We deduce all necessary information from the orthogonal 
polynomials. Consider two families of polynomials: 



Pn{x] t) = w„ ^a;; AVt, BVi, 
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with Wn defined by (|1.15|) . and Q„ that are essentiahy 

Wn (^y, AVi, BVi, ^{x + y^x"^ - 1), - V - 1] 



muhiphed by 



To avoid singularities that may arise when x is in the discrete part of Us, formally 
we define Qn for all x by Q-i = 0, Qo = 1 and by the three step recurrence: 

(3.12) 2yQn{y\x,t,s) 

= AnQn+l{y]X,t,s) + (a + 1/a - An - Cn)Qn{y;X,t,s) + CnQn-liV, X,t, s) 

with An, Cn defined by ((LT6)) and pTf| , 

(l-abcdg"-i)(l~a6g") 
^ ^ " a(l-abcdq^"-^){l-abcdq^")' 

(3.14) 

- _ a(l - q")(l - acg"-i)(l - adq"-i)(l - bcq"-^){l - bdq"-^){l ~ cdq"-^) 
" " (1 - a6cdg2n-2)(i abcdq^''-^) ' 



and with a = A^/i, b = B\/t, c ~ y/j{x + ^/x"^ — 1), and d ~ ^/j{x — \/x'^ — 1). 

Polynomials {pn} are orthogonal with respect to TTt{dx) and polynomials {Qn} 
are orthogonal with respect to Ps^t{x,dy) when the parameters take appropriate 
values. 

The connection coefficients are due to Askey and Wilson [3]. 
Lemma 3.5. For n>l, 

n 

(3.15) Qn{y;x,t,s) = ^bn.k{x,s) {pk{y;t) - pk{x;s)) , 

fc=i 

where bn.kix, s) do not depend on t and bn^n{x, s) 7^ does not depend on x. 
Proof. An adaptation of 4, Formula (6.3)] with j3 = b shows that 

n 

(3.16) Qn{y; X, i, s) = ^ bn,k{x, s)pk{y; t) 

k=0 

with the connection coefficients 

bnAx,s) ^ {-ifq^^^-^^'Ha,, aS)n ^^T/T^ ^^'^ 

[q,q '-abcd,ad,a'^)k 

( ab^Sq'^^'^^^ ,acq^ , adq^ 
""^^^l^ q^''abcd,a^q\a8q^ ' 

which do not depend on t, as 

ab-16 = ABCD, abed = ABs, 



a'y = A^/s{x + \/ x^ — 1), ac = AC, 



a5 = Ay/s{x ~ \/x^ - I), ad = AD. 
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By ([221, 

On,n[X,S)-[ I) q {q,qn-lABs,AC,AD)n 

iq,q"-^ABCD)„ ^ 
{q,q^^-^ABs,AC,AD)n ^ 

does not depend on x. 

We now observe that p.l2p gives (5i(a;; x, s, s) = 0, as a+l/a— Aq— Co = 2x when 
t — s, and Q2(x] x, s, s) — 0, as (71(0, b, c, d) = when t = s. So Qn{x; x, s,s) — 
for n > 1, and (|3.16p imphes 

n 

'^bn,k{x,s)pk{x;s) = 0. 

fe=0 

Subtracting this identity from p. 161) we get p. 151) . □ 

We also need the following generalization of the projection formula [50] . 

Proposition 3.6. Suppose that A, B, C, D satisfy the assumptions in Theorem ] 1 . 1\ 
For X G Us, 



(3.17) J pn{y;t)Ps,tix,dy) = pnix;s) . 

Proof. Since x G Us, from Proposition 13. 3[ measures Ps^t{x,dy) are well defined. 
The formula holds true for n = 0. Suppose it holds true for some n > 0. By 
induction assumption and orthogonality of polynomials {Qn}, 

= / Qn+i{y;x,t,s)Ps,t{x,dy) 

= bn+i,n+iix,s) J {p„+i{y;t) - pn+i{x;s)) Ps,t{x,dy) 

= bn+i^n+i{x,s) (^j pn+i{y;t)Ps^t{x,dy) - pn+i{x;s)^ . 

□ 

Proof of Proposition \3.Ji\ This proof follows the scheme of the proof of p'. Proposi- 
tion 2.5] To prove p. 101) . let im{V) — J^Ps^t{x, V)ns{dx) and note that by orthog- 
onality, J^pn{x; s)Trs{dx) — for all n > 1. Then from (I3.17p . 

Pn{y;t)fj,{dy) ^ i Pn{y,t)Ps,t{x,dy) ] ns{dx) = / Pnix; s)ns{dx) = 0. 



Since JjgPn{y',t)Trt{dy) = 0, this shows that all moments of nijiy) and TTt{dy) are 
the same. By the uniqueness of the moment problem for compactly supported 
measures, n{dy) — TTt(dy), as claimed. 

To prove (I3.11|) . we first note that for x G Us, Ps,t{x, Ut) = 1; this can be seen 
by analyzing the locations of atoms, which arise either from the values of A\/t or 
B^/t > 1 or from x being one of the atoms of Us- (Alternatively, use p.lOp .) 
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Fix X eUs and let ^(F) = ^^^Pt,u{v,V)Ps,t{x,dy). Then, by ([336]) for n > 1 
and p.l7p used twice, 



Thus the moments of /i and Ps^u{x, dz) are equal which by the method of moments 
ends the proof. 



3.3. Proofs of Theorem [T7l] and Theorem 11.21 

Proof of Theorem ] 1.2[ If A,B,C,D satisfy the assumptions in Theorem ll.li by 
Proposition 13.41 there exists a Markov process (Yt) with orthogonal polynomials 
{pn{x;t)}. From p.l7p we see that {p„(a;;i)} are also martingale polynomials 
for (Vt). With Zt defined by (|2.22p . polynomials rn{x;t) inherit the orthogonal 
martingale property, as r„ {Zt;t) = p„ {Yt ■,t). □ 

Proof of Theorem From Proposition 13.41 we already have the Markov process 
(It) with Askey- Wilson transition probabilities. Thus the mean and covariance are 
(|2.14l) and (|2.16p . Formulas (|2.22p and (I2.27|) will therefore give us the process 
(Xt)tgj with the correct covariance. 

It remains to verify asserted properties of conditional moments. Again transfor- 
mation (j2.27p wiU imply ()1.3|) . provided {Zt) satisfies (|2.24|) . For the proof of the 
latter we use orthogonal martingale polynomials (|1.20|) . Our proof is closely related 
to [SI Theorem 2.3]. We begin by writing the three step recurrence as 

xr^i {x ;t) = {ant + I3n)rn+i{x;t) + {jnt + Sn)r{x;t) + {ent + (pn)rn-i{x;t) 

which amounts to decomposing the Jacobi matrix of {r„(x; t)} as tx + y. From 
(|1.15|) with a = A^/t, b = B^/i, c = C/Vi and d = D/Vi we read out the 
coefficients: 



a„ = -ABq"l3n , 

(1 - ABCDq''-^) (1 - ACg") (1 - ^L>g") 
^" " Ay^T~^{l - ABCDq^^) (1 - ABCDq^'^-^) ' 

A {I- BCq"^-^) (1 - BDq''-^) (1 - g") 
^" " yr^(l - ABCDq^"'^) (1 - ABCDq^'^-^) ' 

7„ = A/ \/l-q - a„ - e„ , 




&n,fc(a;, s) {pk{z; u) - pk{x; s)) Pt,u{y, dz)Ps,t{x, dy) 



I y2 ^n-kix, s) {pk{y; t) - pk{x; s)) Ps,t{x, dy) 
u*k=i 




□ 



-5, 



1 



13,1 - ■ 



n 



A^/T^ 
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A calculation verifies the g-conimutation equation [x, y]q — I for the two compo- 
nents of the Jacobi matrix. In terms of the coefficients this amounts to 

(3.18) q;„/3„_i = qan-if3„ , 

(3.19) /3 

(3.20) 7„(5„ + /3„e,i+i + a,i-i0„ = g(7„(5„ + /?„_ie„ + Q!„0„+i) + 1 , 

(3.21) (5„e„ + 7„ _!(/)„ = g(5,i-ie„ + 7„0„) , 

(3.22) = gen+10n . 

For a similar calculation see j30[ Section 4.2]. A more general q-commutation 
equation [x, y], = I + 6'x + lyy + rx^ + cry^ appears in [51 (2.22)-(2.26)]. 

For compactly supported measures, conditional moments can be now read out 
from the properties of the Jacobi matrices: formula (|1.2p follows from Jt = ix + y, 
and formula (II. 3p follows from [x, y\q = I. This can be seen from the proof of O 
Lemma 3.4] but for reader's convenience we include some details. 

Denote by r{x;t) = [ro{x;t),ri{x;t), . . .]. Then the three step recurrence is 
xr(x] t) = r{x; t)3t and the martingale polynomial property from Theorem 11.21 savs 
that E{r{Xu; u)\J-t) — ^{Xt;t). (Here we take all operations componentwise.) 

To verify (|1.2|) for compactly supported measures it suffices to verify that 



u — t t — s 
E{Xtr{Xu;u)\J^s) - E(X,r(X„; u)| J",) + E(X„r(X„; ii)| J",). 

U — 3 U — S 

Using martingale property and the three step recurrence, this is equivalent to 
r{Xs; s)3t = r(X,; s)i^Js + H^-) ^^i^h holds true as J* = ^J, + for 
linear expressions in t. 

To verify (|2.24p we write it as 

(3.23) 

2|^^ J ^ {u-t){u-qt)Z^ ^ iq+l)it-.s)iu-t)ZuZ, ^ {t - s){t - qs)Zl 
* [u—s)(u — qs) [u—s)[u — qs) [u — s){u — qs) 

{t-s)[u-t) 



qs 



For compactly supported laws, it suffices therefore to verify that 

(3.24) E(z2r(X„; u)| J-,) = ~ ^j^" ~ E(ZM^.; ")l-^.) 

[u — s)(u ~ qs) 

, (g + l)(t-s)(u-i) {t~s){t-qs 
-¥.{ZuZsr{Xu]u)\Fs) 



(u — s){u — qs) ' {u — s){u — qs) 

^ (t-.)(u-t) , 
u — qs 

Again, we can write this using the Jacobi matrices and martingale property as 

3.25 r(X,;s)Jt ^r[Xs;s)[- -J, H r- ^ — J,.J„ 

\[u — s)[u — qs) [u — s)[u — qs) 

{t-s){t-qs) , ~ ^^^^ ^ %(X ■ s) 



{u — s){u — qs) "/ u — qs 
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A calculation shows that for J* = ix + y, the g-commutation equation [x, y]q = I 
is equivalent to 

CiOa^ t2 _ - - g^) t2 , (g + - ^)(" - T T . ft - ■'^)(^ - g'^) t2 

\u — s)\u — qs) [u — s)(u — qs) [u — s)(u — qs) 

{t-s){u-t) 
u — qs 

so ([5:^ holds. 

Uniqueness of {Xt) follows from the fact that by [HI Theorem 4.1] each such 
process has orthogonal martingale polynomials; from martingale property p.l7p 
all joint moments are determined uniquely, and correspond to finite dimensional 
distributions with compactly supported marginals. 

□ 

We remark that the following version of Propositions 12.61 and 14.51 would shorten 
the proof of Theorem ll.il 

Conjecture 3.1. Let {Yt)tei be the Markov process from Proposition 13.41 Then 
for any s < t < u from I{A, B, C, D), the conditional distribution of Yt given J^s.u 
is 



(3.27) iy[y; 

(Here, x^Ys + VW^, z^Y^ + ^/W^.) 

4. Purely discrete case 

Assumption \1A\ arises from the positivity condition for the Askey- Wilson re- 
currence for which it is difficult to give general explicit conditions. The following 
result exhibits additional quadratic harnesses when condition (|1.4p is not satisfied. 

Proposition 4.1. Suppose q, A, B,C, D > and ABCD < 1. Suppose that there 
are exactly two numbers among the four products AC, AD, BC, BD that are larger 
than one, and that the smaller of the two, say, AD is of the form 1 /q^ for some 
integer N > 0. If Aq^ > 1 then there exists a Markov process (^t)fe(o.oo) with 
discrete univariate distributions supported on N + 1 points such that (jl.ip (|1.2p 
hold, and (jl.3l) holds with parameters ri,9,a,T,j given by (|1.5p through (11.91) . 



After re-labeling the parameters, without loss of generality we will assume that 
< A < B, < C < D, AC < 1, BC < 1, AD = 1/q^, so that BD > 1/q^. 

Remark 4.1. Discrete Askey- Wilson distributions iy(dx;Ay/i,By/i,C/Vt,D/\/t) 
may be defined also when ABCD > 1, a case which we do not consider in this 
paper. 

4.1. Discrete Askey- Wilson distribution. The discrete Askey- Wilson distri- 
bution i'n{dx;a,b,c,d) arises in several situations, including the case described in 
Lemma [3.1f iii). This distribution was studied in detail by Askey and Wilson [3], 
and summarized in 4 . 

Here we consider parameters a,b,c,d > and < q < 1 such that ad = l/q^ 
and 

q^a > 1, q^ /{be) > 1, ac < 1, q^ a/b > 1, q^a/c > 1, q^ab > 1. 
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Note that this imphes abed < 1 and 

ad = > 1, ac < 1, be < I, bd < 1, ed < 1, ab > 1/(J^, 

so from Lemma IS-lT iii). Askey- Wilson law i'{dx;a,b,e,d) — i'{dx;a,b,e,l/{aq^)) 
is well defined, depends on parameters a, b, e, q, N only, and is supported on + 1 
points: 

{Xk = {q'^a + q~''a-^) /2 : fc = 0, . . . , TV} , 
According to [3], the Askey- Wilson law assigns to x^ the probability Pk,Nio-, b, e) — 
Pk{a,b,e, l/{q"a)) (recall p.GI) '). The formula simplifies to 



(4.1) 

Pk,N{a, b, c) 



^),_,(«^ «c),(l-,-«^),^ 



Here 



{q'^a^)N+l 
denotes the g-binomial coefficient. 



0,. 



{q)k(q)N-k 

We remark that if A is a random variable distributed according to v{dx; a, 6, c, 1/ {aq^)) 
then E(A) and Var(A) are given by formulas (1^ and (P?71) with d = l/(q"a) re- 
spectively. This can be seen by a discrete version of the calculations from the proof 
of Proposition l2.2l alternatively, one can use the fact that the first two Askey- Wilson 
polynomials wi{X) and W2{X) integrate to zero. 

The discrete version of Proposition 12.31 savs that with d = l/{maq^), 

(4.2) i'{U;a,b,e'm,md) 



iy'{U; a, b, ra{x + \J x"^ — 1), m{x — x'^ — \)^v(dx\ ma, mb, e, d), 
and takes the following form. 

Lemma 4.2. For any m E (0, 1) and any j = 0,1, . . . , N 



N 



(4.3) 



PjM{a, b, me) = ^ P],k{a, b, q^m'^a)pkM{ma, mb, e) 

k=3 



Proof. Expanding the right hand side of (|4.3p we have 



N 

(4-4) E 

k=j 



b ' q^-^m2j^_j [ab, q m^a^)j {1 — q^^ a^)q 



{qla^) 



k+1 



y q^m? ab J ^ 



(—q'^m'^aby 



"A^I V & ' c ) N~k {rn^ab, mae)k {I — q^^ m^ a^)q~ 



(g'=m2a2)jv+i {:;^^ ^ {-mbeY 

(4^-^). {ab,mae),{l~aW^)q^ 



{<i'^^)N+i{^)^{-mcby 



N 

E 

k=j 



k~ j 



ma k+1 2 

q^ a 



N-k 



{m^)k-j{mq'ac-)k-j{l - q^''m^a^)q 



(fc-j)(fc-3 + l) 



k-3 
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Here we used identities ()2.ip and (|2.2p . The first one for: (i) a = q^^^a/b, 
M = k-j and L = N ~ k] (ii) a = q^m^a^, M = j, L = N-j + l; (iii) a = g-^a^, 
M = k + 1, L^N — k. The second one for: (i) a — m?ah and M — k; (ii) a = in? 
and M — k ~ i . 

We transform the sum in (j4.4p introducing K — k j , L — N — j , a — mq^a, 
/3 = and 7 = c. Then by dUl]) we get 



E 



,-+^f)^_^(a/3),_, (a7)K (1 - ,-a^),^ 



-/37)^ 



= 1 . 



Now the resuh fohows since the first part of the expression (|4.4p is the desired 
probabihty mass function. 

□ 

4.2. Markov processes with discrete Askey- Wilson laws. We now choose the 
parameters as in Proposition liH 0<yl<B, 0<C<i:' = l/{Aq^), ABCD < 1, 
EC < 1, and choose the time interval / = {C{q'^A)-^, (AB)-^) from pH^ . For 

any t £ I, define the discrete distribution TTt{dx) — J2k=o'^tiyki^))^yk{t){dx) by 
choosing the support from (13. 9p with weights 



(4.5) nt{ykit))^Pk,N[Asi,Bt^,Ct-^^^ . 

Clearly, the support of tt^ is 11^ = {yo{s), yi(s), . . . , 2/jv(s)}. 
Also for any s,t G I, s < t and for any k G {0, 1, . . . , A^}, define the discrete 
Askey- Wilson distribution Ps^t{yk{s), dx) = J2j=Q Ps,t,yk{s)^yj(t){dy) by 

(4.6) P,^t^y^^s){ydt))=PjM{Ati,Bti,q''Ast-^'^ . 

Thus Ps,t{x,dy) is defined only for x from the support of tt^. Next, we give the 
discrete version of Proposition 12.41 



Proposition 4.3. The family of distributions {nt, Ps.t{x,dy)), s,t £ I , s < t, 
k e {0, 1, . . . , A'^} defines a Markov process {Yt)t£i with trajectories contained in 
the set of functions {{yk{t))tei, fc = 0, 1, . . . , N}. 

Proof. We need to check the Chapman-Kolmogorov conditions. Note that for any 
s < t and any k the support of the measure Ps,t,yk{s) is a subset of the support Ut 
of the measures nt. First we check 

AT 

(4.7) MvAt)) = E Ps,t,k{yj{t)) T^s{yk{s)) 

k=j 

which can be written as 

N 

PjM {At^,Bt^,Ct-^^ = ^ pj^k {m^- ,Bt^ ,q^ Ast-^^ pkM {As^,Bs^,Cs-^ 

k=3 

Now (|i?7)) follows from with a = Ati , b = Bt^ , c = Cs~i and m = (s/i) - 
Similarly, the condition 

i 

(4.8) Ps,u,u^:{s){Ut{u)) = ^ Pt,u,u,{t)(ui{'U-)) Ps,t,Uk(s){Uj{t)) 

J=k 



ASKEY-WILSON POLYNOMIALS, QUADRATIC HARNESSES AND MARTINGALES 27 



assumes the form 

i 

= pij{Au^ ,Bu^ ,q' Atu''^)pj^k{At^ ,Bt^ Ast'^) . 

j=k 

Therefore (|4.8p follows from (|4.3p with {j,k,N) {i,j,k), a = Au^ , b = Bu^ , 

c = q'^Ast^^ and m ~ {t/u)^ . □ 

Let {Yt)tei be a Markov process defined by the above Markov family {tti, Ps,t.yk(s))- 
Note that at the end-points of /, i^c/(g"A) is degenerate at ^^^^"^ i and yi/(^B) 
is degenerate at ^+^^^ (compare Proposition 15 . 1 f i) ) . 

2{AB) 2 

Proposition 4.4. For the process {Yt)tei defined above 

^ (1 - q^)A{Bt + C)-{1- q^AH)il - BC) 
^ " 2At-2iBC-q^) 

Cov(i;,y() 

_ (1 - g)(l - - AC)(1 - BC)(g"A - B){q^As - C)(l - ABt) 

^ Aq{st)i{BC - q^Y{BC - q^'^) ' 

^ (1 - AC){1 ~ BC){1 - ^ g)(l - ABt){C - q'^At) 

^ " iqAt{BC ~ q^)^BC ^ q^-^) 

Expressions for conditional expectations and conditional variances are exactly the 
same as in the absolutely continuous case. 

Proof. The result follows from the fact that the marginal and conditional dis- 
tributions of the process {Yt)tei are finite Askey- Wilson. Therefore one can ap- 
ply the formulas (|2.6I) and (|2.7p . The covariance is derived through conditioning 

EiYsYt)^E{YsEiYt\Ts)). □ 

Note that the above formulas are identical to those in the case with densities by 
plugging there D = q^^ A^^. 

Since we are interested in the harness properties, we want to find the conditional 
distributions of the process with conditioning with respect to the past and the 
future, jointly. The following result says that the conditional distribution of Yt 
given the (admissible) values Yg = x,Yu = z is a discrete Askey- Wilson distribution 

1/ ^dy; ^{z + Vz^-l), \l^{z - y/z^ - 1), -^-{x + y/x^ - 1), ^-{x- x^ - 1)^ , 

compare Proposition 12.61 (Using notation (|4.ip . this formula takes a much more 
concise form.) 

Proposition 4.5. Let {Yt)t£i be the Markov process defined by (ttj, Ps,t,yk{s)) given 
by (j4.5p and (|4.6p with parameters A, B, C, q, N. Then for any s,t,u € I such that 
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s < t < u the conditional distribution of Yt given Fs,u is defined by the finite 
Askey-Wilson distribution 



P{Yt = y,{t)\Ys - yk{s),Yu = y^{u)) = PJ-^M-^ q'At^, 



t2 q'^As 



The expressions for the first two conditional moments are the same as in the abso- 
lutely continuous case. 

Proof. Due to the Markov property of the process {Yt)t(zi, to determine the condi- 
tional distribution of Yt given Fs.u it suffices to find P(Yt = %(i)|^s = yfe(s), F« = 
yi{u)) for any i, j, k £ {0,1,..., N}, such that i < j < k. Also the Markov property 
implies that this probability can be expressed in terms of conditional probabilities 
with respect to the past as 

p{j\k,i) - P{Yt = yj{t)\Ys = yk{s),Yu = y^{u)) 

^ P{Y^ ^ y,{u)\Yt ^ y,{t))P{Yt = y,{t)\Y, = y^js)) 
P{Yu ^ y^{u)\Ys ^ ykis)) 

_ Pij{Au^ , Bu^ ,q^ Atu^i) pj^iAt^ , Bt^ ,q''Ast^i) 
Pi^k{Au^ , Bui , q''Asu~i) 
Expanding the expression for the probability mass functions according to (|4.ip 
we get 



qu/t\ 

{ B ' {ABu,q^AH), {l-q^'A^u)q— 



{q'A^u) 



( q \ {-qiABty 

\q^ABtJ. 



( q' + ^A qt/s\ 



{ABt, q''A\s)^ (1 - q^^AH)q'^ 



{q^ABs), 



-q'^ABsy 



f q"-^^ A qu/s 
\ B ' q''-' 



k^, {ABu, q'^A^s), (1 - q^''A^u)q 



{q'A^u)k^ 



1^'ABt 



{-qf'ABsf 



This can be reduced in several steps. The g-binomial symbols reduce, as in the 
^ * ' Then we apply (12. ip in the following situations: 



classical {q—1) case to 

(i) a = 2!^^ M=j-i,L = k~j, (ii) a = q'A^u, M = j + l,L = k- j, (iii) 
a = q^ A?t, M — i, L = k — i + 1, (iv) a = q'^A'^s, M = i, L = j — i. Also we apply 
([2:21) for (i) a = f , A/ = j - i, (ii) a = ABt, M = j. Thus 

U-i)U-i + i) 



pU\k,i) 



k — i 
j - i 



ipl, q^+'^+'A^u) . {I q'^+^A^s) (1 - q^^AH)q''-^ 



{q^+iAH)k- 



k—i 



which, through comparison with the definition ()4.ip , is easily identified as the dis- 
tribution we sought. 

□ 
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Proof of Proposition \41\ Since formulas ()2.6|) and (|2.7p hold for all Askey- Wilson 
distributions, from Proposition 14.51 we see that the conditional moments and vari- 
ances in the discrete case are also given by formulas from Proposition l2.5l Therefore 
the transformed process 

(1 + ABt)2T(t)iYT(t) -iA + B)t - (C + -4^) , 

^(1 - - AC){1 q-^){l BC){1 - ^) 

is a quadratic harness on J with 9, r], r, a, 7 defined as in the general case with 
D = q-^A-^. □ 

5. Concluding observations 
This section contains additional observations that may merit further study. 

5.1. Bridge property. The following proposition lists some combinations of pa- 
rameters that create a "quadratic harness bridge" either between two-point masses, 
or degenerated laws. 

Proposition 5.1. Let {Zt)tei be the Markov process from Theorem ] Assume 
that AB ^ so that (|1.19p defines a finite interval I ~ [So,Si], and extend Zt to 
the end-points of I by L2- continuity. 

(i) If AB > then ~ (1/^+ l/i3)/\/l — q is deterministic; similarly, if 
CD>0 then Zs, = {C + D)/y/T^ . 

(ii) If q < and CD < then Zs-^ takes only two-values. Similarly, if q < 
and AB < then Zs^ is a two-valued random variable. 

(iii) If CD < and q > then Zq is purely discrete with the following law: 

(^^^ V (7 iAD,BDUiAC,BC)k ^ ^ ^ 

^'■'^ '^'r=VT=q) = iD/C,ABCDUiq,qC/D),'^^'^'^ 
q'^D \ _ {AC,BCU{AD,BD)k 



(5-2) P\Z,-^==j - ^c/D,ABCDU{q,qD/C),''''''-^- 

Proof. We can derive the first two statements from moments, which are easier to 
compute for (Yj) instead of (Zt). In the first case, Var(yt) = at the endpoints, 
see (|2.15l) : in the second case i?(ui|(yt)) = at the end-points. Alternatively, one 
can compute the limit of the Askey- Wilson law as in the proof of part (iii). 

For part (iii), without loss of generality assume |^| < \B\ and \C\ < \D\. Then 

the discrete part of Zs has atoms at | ^i^q (q''^ + c^) ■ J — q^-'C^ > s| and 

The probabilities can be computed from p.6p with c ~ Ay^, d — B^ and either 
a = C/y/s, b — D/y/s for (|5.ip or a = D/y/s, b — C / y/s for (|5.2I) . and converge to 
(|5.ip and ()5.2p respectively. To see that the limit distribution is indeed discrete, 
we note that 



00 



fe=0 



{AD,BD)a, /AC,BC \ {AC,BCU (AD,BD 

{DICABCDU V 9^/^ ''^r{ClD,ABCDU V <lD/C ' 
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Here we use hypergeometric function notation 



Oi, 02, . . . ,ar+i \ _ir-^ (ai,a2,. . ■ ,ar)k k 



/ "1 , "2 ; ■ • ■ 7 «r+l \ _ ST^ 



z 



A = 



has E(XtXs) = min{s,t}. The easiest way to see this is to note 



The identity that gives the final equahty is [THl (12.2.21)] used with a — AC, b — 
BC,c^qC/D. □ 

5.2. Transformations that preserve quadratic harness property. The ba- 
sic idea behind the transformation p.27p is that if covariance K{ZtZs) ~ cq + 
ci min{i, s} + C2 max{t, s} + csts factors as (s — a)(l — t/3) for s < t with a/3 < 1 
then it can be transformed into min{t, s} by a deterministic time change and scal- 
ing. 

This transformation is based on the following group action: if ^ = ^ ^ G 

Gi2(IR) is invertible, then A acts on stochastic processes X — (Xt) by ^(X) :— Y 
with Yt — [ct + d)XTj^{t), where TA{t) = {at + b)/{ct -\- d) is the associated Mobius 
transformation. It is easy to check that this is a (right) group action: A{B{1C)) = 
{B X A)(X). 

If E(ZtZs) = (s - a)(l - tp) for s < i and aP < 1 then X = ^-^(Z) with 
1 -a] 
-P 1 J 

that Ta is increasing for a/3 < 1, and by group property Z = ^(X). So with s < t, 
nZsZt) = (l-s/?)(l-t/3)E {XT,(s)XT,(t)) = {l-sp){l-tP)TA{s) = (5-a)(l-t/3). 

It is clear that, at least locally, this group action preserves properties of linearity 
of regression and of quadratic conditional variance. In fact, one can verify that the 
general form of the covariance "E^^Xt, Xg) = co -I- ci min{i, s} -I- C2 max{t, s} -I- csis 
is also preserved, and since this covariance corresponds to (|1.2I) . the latter is also 
preserved by the group action. 
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